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A conspicuous tangential alignment of galaxies in a STIS 
Parallel Shear Survey field: A new dark-lens candidate ? * 
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Abstract. We report the serendipituous discovery of a conspicuous alignment of galaxies in a field obtained through the STIS 
Parallel Shear Survey. This project collects randomly distributed 50" x 50" fields to investigate the cosmic shear effect on this 
scale. Analyzing the parallel observations having the Seyfert galaxy NGC625 as primary target, we recognized over the whole 
field of view a strong apparent tangential alignment of galaxy ellipticities towards the image center. The field shows several 
arclet-like features typical for images of massive galaxy clusters, but no obvious over-density of bright foreground galaxies. We 
also find a multiple image candidate. On the basis of the possible strong and weak lensing effect within the data, we discuss 
whether this could be compatible with a massive halo with no clear optical counterpart. 
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1. Introduction 

The STIS CCD camera on-board Hubble Space Telescope has 
been providing an unique view of faint galaxies. The charac- 
teristics of this instrument which offers good depth, excellent 
resolution and a good sampling of the HST PSF, together with 
the possibility of obtaining data using the parallel observing 
mode, makes it a suitable instrument for studies such as the 
STIS Parallel Shear Survey (hereafter refered as SPSS), whose 
main goal is the study of cosmic shear on scales below one 
arcminute. Data for this program started to be collected on 
September 24th, 2000, within the scope of an HST GO cycle 9 
parallel proposal (8562 and 9248, P.I.: P. Schneider) and up to 
date, more than 500 different pointings have been obtained. The 
use of STIS CCD imaging in the CLEAR filter mode allows 
to combine good spatial resolution (originally 0.05" per pixel) 
with very deep exposures since the CCD is un-filtered and sen- 
sitive to wavelengths between 2500 and 1 1000 A(Pirzkal et al. 
2001). We report in this paper the discovery of a conspicuous 
alignment of galaxies in a deep exposure obtained with that 



program. This field contains several high-ellipticity galaxies 
which seem oriented around a point near the center of the field 
giving a striking visual impression that caught our attention. 
Also, there are several fainter galaxies in this field which also 
contribute to the visual impression of an alignment of galax- 
ies. We used the aperture mass technique and an analysis of 
the most elliptical objects in the field to try to quantify the sig- 
nificance of this effect and test the visual impression. Section 
2 describes our data and the reduction procedure. The core of 
section 3 presents the aperture mass analysis, discussing the 
significance of the peak of mass obtained and the errors on its 
position. Section 4 presents the most elongated objects in the 
field, and discusses their orientation with respect to the peak 
found with the aperture mass technique. We discuss also in this 
part the similar characteristics of 2 of these candidate arclets 
that could be compatible with being multiply imaged. Finally, 
section 5 summarizes and discusses the indications presented 
in the paper which, taken together, could make a case for the 
presence of a possible underluminous lens in this field. 
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* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Space Institute (STScI) 
which is operated by AURA Inc., under NASA contract NAS 5-26555 



These data were obtained as part of the large HST GO cycle 
9 parallel proposal ai med at the study of cosmic shear ( see 
Pirzkal et al. 2001 and http://www.stecf.org/projects/sheai for 
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more details on the STIS data characteristics). The field con- 
sidered here was, coincidentally, the first one observed as part 
of that program. The associated image is composed of 18 in- 
dividual 400s exposures with a field of view of 50" x 50", 
taken in the CLEAR filter mode. The exposures were individ- 
ually cleaned for hot pixels and cosmic rays using the Eye al- 
gorithm in SExtractor (Bertin & Arnouts 1996)[], which cre- 
ates masks of the pixels that are not used during the combi- 
nation, eliminating 99% of the spurious pixels. The shifts be- 
tween the different exposures were computed using the method 
described in Pirzkal et al. (2001), which can align these im- 
ages to an accuracy of 1/10 of a pixel. Finally, the individual 
exposures were drizzled to a pixel size of 0.025" (Fruchter 
and Hook 2002) and combined by median averaging with the 
IRAF Imcombine procedure, achieving a total exposure time of 
7200s. The use of this reduction procedure, basically similar to 
the one described in Pirzkal et al. (2001), is known to preserve 
the original PSF shape and size (Hammerle et al. 2001) which 
is extremely important for an accurate analysis of the shear. 
The coordinates of the field are R.A: 01 h 35'06.7" and Dec: 
-41°21'20.18" (J2000), about 7' north of the local Seyfert 
galaxy NGC625, which was the primary target observed with 
the WFPC2. Throughout this paper we will refer to this field 
as the Slensl field. The coadded image is available in fits for- 



mat at frttp ://w w w. stecf.org/proj ec ts/shear/slens/slens . fits . The 
SExtractor parameter file used is also available there in order 
to enable others to reproduce the SExtractor catalog. 



rection (see also Hoekstra et al. 1998 for a similar procedure 
with WFPC2 data). The final correction was performed using 
the PSF information from the starfields as described in Erben 
et al. (2001). We removed from the catalog the objects with an 
axis ratio larger than two; those will be considered in the next 
section as candidate arclets. Their alignment towards the direc- 
tion defined by the A/ ap statistics will provide an independent 
test of the shear signal. In the end, our catalog contained 54 
objects, with fully corrected ellipticities, which were used for 
the subsequent analysis. 




3. M ap statistical analysis 



Over the whole field of view of our 50" x 50" STIS image, 
we recognize a strong apparent tangential alignment of galax- 
ies towards the center (see Figs, [j] and ||). Although the small 
50" x 50" field does not allow a quantitative weak lensing anal- 
ysis to verify or to falsify the hypothesis of a massive dark mat- 
ter halo in the image center, we use the M ap statistic (Schneider 
1996) to test whether we can identify a significant and sta- 
ble peak. The details of the creation of an object catalog are 
the same as for the rest of the SPSS project. They are de- 
scribed and justified in depth in Hammerle et al. (2001). In 
short, we performed an analysis of the data using two differ- 
ent software packages. Objects were detected with SExtractor 
and a modified version of the IMCAT package (Kaiser, Squires 
& Broadhurst 1995 hereafter referred to as KSB; Erben et al. 
2001). Finally, we only retained objects detected with both 
programs, assigning them the photometric information from 
SExtractor and the shape parameters from IMCAT. There are 
not enough foreground stars in a STIS image to allow directly 
a PSF correction in the KSB style from the data itself, which 
is the usual procedure with wide-field ground-based observa- 
tions. On the other hand, our analysis of STIS data, especially 
starfields obtained over the same period of time as our obser- 
vations, showed that the PSF is sufficiently stable (Hammerle 
et al. 2001) so that these starfields provide a good estimate of 
the quantities used for the KSB anisotropy and smearing cor- 



Fig. 1. The figure shows contour maps of the S/N of M ap val- 
ues (eq. ^) superimposed on the field. The S/N ratio was cal- 
culated on a 20x20 grid of our image. The filter scale 9 is 25". 
White contours correspond to S/N = 1.0, 1.5, 2.0 and black 
contours to S/N = 2.1, 2.2, 2.3. See the text for a discussion 
on the robustness of this result. 



We use standard lensing notations in this paper. For a 
broader introduction to this topic, the reader should refer to 
Bartelmann & Schneider (2001). We used the M ap statistic 
(Schneider 1996) for the weak lensing analysis. The family 
of M ap statistics uses the fact that a filtered integral over 
the convergence k can be converted into a filtered integral 
over the observable tangential shear j t [see equations (9) and 
(15) in Schneider 1996], if the filter function J7(|i?|) satisfies 
J AO' 6' U{6') — but is arbitrary otherwise. It is straightfor- 
ward to construct an unbiased estimate M ap for the integral by 
a discrete sum over observed image ellipticities e t and consid- 
ering the coordinates origin being at the center: 



7T<9 



1 Eye an d SExtractor are part of th e Terapix software suite and can 
be found at ittp://terapix.iap.fr/soft/| 



with Q{x) = q{x/9)/9 2 and q(p) = (6/tt) p 2 (1 - p 2 ). 



(1) 
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Moreover, as it is a scalar quantity, expressions for the vari- 
ance and hence the S/N for a derived M ap value are easily 
calculated. For the S/N one obtains (see Schneider 1996): 

s= V2 ZM0i)Q(\0i\) 

where a t is the intrinsic dispersion of image ellipticities. It is 
estimated by the dispersion of the galaxies entering the calcula- 
tion of S [of = N^ 1 J2i( e ii + e 2i)l- Hence the M ap statistics 
are the ideal tool to evaluate the significance and robustness of 
a mass concentration. In our case, with only a 50" x 50" field 
and 54 objects in hand, the possibilities for this kind of analysis 
are limited. We are confined to small radii for our filter scale 
9 (about 25") for which the M ap statistic is extremely noisy. 
Moreover, applying ([j]) to a grid of points on our data will not 
give an unbiased estimate of M ap since the data do not cover a 
complete circle with radius 8 at most positions. However, even 
if M ap can no longer be interpreted as an unbiased estimate of 
the filtered density profile, it yields a measure of the tangen- 
tial alignment around a point and as such is a useful statistical 
quantity. 

The application of this approach to our data gives the re- 
sults shown in Fig. |J There, we show a significance map for 
one filter radius (6 = 25"). We could slightly vary the size of 
the filter, but we are limited by the size of the field and the 
number of objects in it. Changing the filter size down to 20" 
or up to 30" does not modify the results. We recover a peak 
located at the image pixel coordinates x — 1066; y — 1297 
(here and in the following, 'pixels' denote the subsampled pix- 
els of size 0.025" and with an origin for the coordinate system 
at the lower left corner of the image), whose position approx- 
imately coincides with the position of the brightest galaxy in 
the field. With a S/N s» 2.42 for this smoothing, the peak 
is marginally significant. Here, the significance was estimated 
with (||). We found that these estimates are in very good agree- 
ment with those obtained from randomizing the orientations of 
galaxy ellipticities. To judge whether the peak is a robust sta- 
tistical signal, we repeated the calculation of the M ap statis- 
tics drawing randomly 70% and 50% of the galaxies out of 
the initial 54 objects catalog. This exercise was repeated 100 
times and the positions of the most significant peaks in the re- 
alizations were considered. The mean and standard deviation 
for the recovered positions of the peak are in the 70% case: 
(x) = 1250; a x = 265; (y) = 1383; a y = 244 and for the 
50% case: (x) = 1234; a x = 300; (y) = 1286; a y = 350. The 
number of "catastrophic outliers", where the most significant 
peak is located at the very border of the frame, is about 12% 
in both cases. The mean amplitude for the peaks is S/N ss 2.4 
(70% case) and S/N « 2.1 (50% case). These values are a bit 
higher than expected from simple Gaussian statistics. Also if 
we choose our samples not randomly but remove the 30% and 
50% most elliptical galaxies from the initial catalog, we recover 
peaks consistent in position with these results (although in the 
30% case, a strong peak appears at the border of the frame). 

Finally we checked the probability of obtaining at least a 
S/N « 2.42 detection in fields with randomly oriented galaxy 
ellipticities. As mentioned above, equations (|IJ) and (||) applied 



to our field do not provide unbiased estimates of Af ap and its 
uncertainty, since our data do not cover a complete circle of 
the smoothing radius at the peak position. In order to estimate 
this probability, we generated 10000 50" x 50" fields with 
54 objects randomly placed, with random galaxy ellipticities 
drawn from the distribution P(ei,£2) oc exp(— |e| 2 /cr 2 ) with 
er c = 0.36 matching the a t of our data. We find that we can 
reproduce a peak of at least S/N w 2.42 with 700 < £ pca k < 
1300; 1000 < ypcak < 1600 in 2.85% of all cases. We con- 
clude that the 54 objects detected in our frame show a fairly 
robust alignment towards the upper middle part of our image. 
As stated above, the lack of knowledge about the shear outside 
the Slensl field does not permit decisive conclusions about the 
existence of a possible massive structure within the field. 



4. Arclet-like features and multiple object 
candidate. 

We selected arclet candidates solely on the basis of their elon- 
gation. The criterion choosen was that the elongation, as de- 
fined in SExtractor (major axis over minor axis), should be 
greater than two. Eleven objects met the requirement, with 
isophotal magnitudes ranging from 22.7 to 27.3 (see Table 
1), which are shown in Fig. 2. A2, which appears to have 
2 components, is detected by SExtractor as two distinct ob- 
jects (A2a and A2b), each with axis ratio larger than two. As 
stated in section 3, these objects were not included in the cata- 
log used to compute the M ap statistic. Eight of those galaxies 
are distributed tangentially to the position of the peak center 
at x = 1066; y = 1297, but three of them (A5, A6 and A8) 
are oriented more radially. The angle between their minor axis 
and the direction towards the Af^ p -peak center is on average 
24 degrees. The probability that the mean of this angle for 11 
objects is smaller than the measured one form randomly ori- 
ented objects is only 0.3%, indicating that those galaxies are 
well oriented tangentially with respect to the peak center that 
was defined by the less elongated objects. This seems difficult 
to explain unless we consider a gravitational effect. 

The high-resolution of STIS CCD imaging permits the 
identification of candidate multiply imaged galaxies, even in 
the absence of any color information as it is the case here, 
just by using morphology and surface brightness. The strongest 
case for a multiply imaged galaxy candidate comes from ob- 
jects Gl and G2 (Fig. 2). These two galaxies show an unusual 
and similar V shape morphology, superimposed on an extended 
and dimmer halo. They have also equal maximum surface 
brightness, 23.29 mag/arcsec 2 for Gl and 23.34 mag/arcsec 2 
for G2. They are located at opposite sides of the M ap peak: 
Gl being at 19.45" from the position of the peak and G2 at 
22.53". Also their tangential alignment relative to the direction 
of the peak is nearly perfect. This symmetry, combined with the 
characteristic morphology and the identical surface brightness 
of these objects, may indicate that they are multiple images of 
the same background galaxy. However, the two images seem 
to show the same parity, with a brighter star-forming region lo- 
cated on the right side, when we would expect an odd parity in 
the case of a double image. 
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Table 1. Characteristics of the objects selected as arclet candi- 
dates. Magnitudes are AB isophotal magnitudes for the STIS 
CCD clear mode filter. Positions are given in image subpixels. 
The elongation is defined by major axis over minus axis. 9 is 
the angle in degrees between the major axis and the X-axis. 



Name Magnitude X Y Elongation 6 



Gl 


23.44 


366 


1697 


2.1 


61 


G2 


23.99 


1870 


1475 


2.5 


-77 


Al 


22.74 


835 


354 


2.0 


-10 


A2a 


23.25 


634 


1961 


3.5 


19 


A2b 


24.12 


679 


1975 


2.6 


17 


A3 


24.24 


287 


942 


2.7 


-46 


A4 


25.76 


1914 


196 


2.4 


58 


A5 


25.86 


1025 


801 


2.5 


56 


A6 


26.11 


1697 


1062 


3.0 


-75 


A7 


26.79 


947 


1454 


2.6 


10 


A8 


26.91 


621 


1503 


3.6 


-62 



Without color information, the question of whether Gl and 
G2 are indeed a double image is impossible to decide. The in- 
dication of the same parity could be false if the bright spot on 
the right in G2 is an unrelated object. However, if they have 
indeed the same parity, a lens interpretation is difficult. Simple 
three-image lens models would predict a third image within 
the field-of-view, unless the putative lens is extremely massive. 
One such example would be a barely critical lens, whose caus- 
tic consists only of a lips singularity (see Schneider, Ehlers & 
Falco 1992); in this case, two equal parity images appear on 
either side of the lens, with the third image near the center. No 
indication for such a central image is found. 

5. Discussion 

The data presented in the previous sections indicate a system- 
atic alignment, and possibly a multiply-imaged object in the 
field. The easiest interpretation of these findings is the exis- 
tence of a mass in the line of sight of Slensl which could pro- 
duce both strong (one candidate multiple image galaxy but with 
a puzzling parity problem) and weak (a coherent tangential 
alignment which can be quantified through the M ap -statistics 
as a significant peak together with several candidate arclets 
distributed tangentially around the same position) gravitational 
lensing effects. However, with our current STIS data, it is not 
possible to proceed further with the analysis and to confirm that 
a massive dark halo is present in this field. The significance of 
the M ap peak is too low to discard the possibility that it is a 
statistical fluctuation, even if we would not expect a better sig- 
nificance given the size of the field and the number of objects 
in it. We therefore considered the result of the M ap -analysis 
only to define a 'center' of alignment, relative to which the 
orientations of the arclets were measured. The mean of these 
orientations is small and is expected to occur randomly with a 
probability of only 0.3%. The arclets show a pattern strongly 
resembling those of lensing clusters. 



If these image features are produced by a gravitational lens, 
it does not seem to be associated with any luminous objects vis- 
ible in the deep exposure obtained with STIS. Only the location 
of the M ap peak could suggest that this mass may somehow be 
associated with the group of bright elliptical galaxies on the 
center-top part of the field. If we suppose that this is the case, 
the mass necessary to produce the candidate multiply imaged 
galaxy like Gl and G2 within a 20" radius would be similar to 
the mass of a rich cluster of galaxies. If the lens is really asso- 
ciated with those bright galaxies located close to the peak, then 
the M/L ratio is at least 2 orders of magnitude above the normal 
values for groups of elliptical galaxies. In any case, if this is a 
lens, the mass we find is unacounted for in terms of light for 
what is normally regarded as the M/L ratio in groups or clus- 
ters of galaxies. We tested also the light profile of the brightest 
elliptical in this group. This galaxy could be at a redshift of 
0.5-0.6 based on similar galaxies in terms of magnitude and 
apparent size for which we have spectroscopic redshifts from 
the CFRS (Lilly et al. 1995). If this galaxy had a cD profile, 
it would strongly suggest the presence of an underlying clus- 
ter, but this galaxy clearly has a deVaucouleurs profile, which 
makes the test inconclusive. 

In conventional models of structure formation, the most 
massive halos detected by aperture mass techniques like the 
one described in section 3 would be easily observed because 
they should have accreted gas which, after subsequent cool- 
ing, would form stars. Nevertheless, it is possible that one 
could find a number of dark lenses, where the dark matter 
would have collapsed in structures massive enough to create 
a strong lens effect, while prohibiting baryonic matter to settle 
in the dark matter halos and initiating the star formation. So 
far, two other candidates have been found (Erben et al. 2000, 
Umetsu & Futamase 2000) but neither of them has been clearly 
confirmed. Also, a class of X-ray Over-Luminous Elliptical 
Galaxies (OLEGs, Vikhlinin et al. 1999) show similar char- 
acteristics of a dark halo, being compatible in X-ray with the 
mass of a poor cluster but presenting a very low optical lumi- 
nosity. Recently White et al. (2001) have suggested that those 
detections could be explained by projection effects of several 
mass sheets along the line of sight. Another explanation could 
be that most of the objects observed in this field are at the same 
redshift and belong to the same gravitational structure, the pre- 
ferred orientation of the galaxies being then just a product of 
their infall into the gravitational potential. 

To confirm the reality and the nature of the M ap peak, wider 
deep exposures centered on the position of the peak are neces- 
sary. At the same time, wider deep multicolor imaging would 
allow us to obtain photometric redshifts that would give an in- 
dication for the redshift of all the objects, including Gl and 
G2, in the field and its surroundings. If Gl and G2 present sim- 
ilar photometric redshifts then, if this is a lens, a third image 
should be detected in those images unless this is a very unex- 
pected configuration. The combination of weak-lensing anal- 
ysis and photometric redshift information can then be used to 
unveil the true nature of the objects in this field. Wittman et al. 
(2001) have recently used the combination of these techniques 
for the discovery of a new cluster of galaxies. The detection 
and confirmation of only a single dark-lens would already have 
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Fig. 2. This figure shows the field with a zoom on the two candidates for multiply-imaged galaxies Gl at the left and G2 at 
the right side of the field and 9 other arclet-like candidates (A1-A8) with isophotal AB magnitudes ranging from 22.7 to 26.9. 
Note that the 2 components of A2 (A2a and A2b) are considered as 2 different objects. The 'X' marks the M ap peak of the 
25" smoothing of Fig. 1. The 2 multiple image candidates are located on each side of the brightest galaxies. They show a very 
characteristic V-shape morphology and have the same surface brightness. 



important consequences for the standard picture of structure 
formation (Trentham et al. 2001). It would also show the im- 
portance of using gravitational lensing as a tool for surveys, 



selecting objects by mass instead of by light as in traditional 
methods. 
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